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Abstract 

The k ■ p method is used to analyze the problem of intervalley T-X^ inter- 
action of conduction band states in (001) lattice-matched III-V semiconductor 
heterostructures. A convenient basis for expansion of the wave function is system- 
atically selected and a multi-band system of equations is derived for the envelope 
functions which is then reduced to a system of three equations for three valleys 
(Fi, Xi, and A3) by using a unitary transformation. Intervalley F-A^ mixing 
is described by short-range potentials localized at heterojunctions. The expres- 
sions for the parameters determining F-A^ mixing strength explicitly contain the 
chemical composition profile of the structure, so mixing is naturally stronger for 
abrupt heterojunctions than for structures with a continuously varying chemical 
composition. It is shown that the direct Fi-Ai interaction of comparable strength 
to F1-A3 interaction exists. This must be taken into account when interpreting 
tunnel and optical experiments since Ai valley is substantially lower in energy 
than A3 valley. 

PACS: 73.20.Dx, 71.25.Cx, 73.40.Kp 



1 Introduction 

Heterostructures with quantum wells and barriers based on GaAs/AlGaAs are popular 
objects for studying the physical processes accompanying resonant tunneling and optical 
transitions. These structures are particularly interesting because in GaAs the minimum 
of the conduction band is situated at the center of the first Brillouin zone (Fi symmetry) 
whereas the minimum of the AlAs conduction band lies near X point where the bands 
having Xi and X3 symmetry are close together, Xi band being the principal one (A3 
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band is 350 meV higher energy). On each cubic axis, there is a pair of X valleys, one from 
Xi band and one from X3 band. Layers with predominantly Ga content act as barriers 
for X electrons while regions with Al predominating act as barriers for F electrons. Under 
certain conditions, resonant tunneling interaction can take place between energetically 
close states of V and X valleys, and this is observed on the current- voltage characteristics 
even of single-barrier structures (see, e.g., Refs. JI], 0]) and also in optical experiments 
The conditions for the existence of this interaction generally involve the presence of 
structural defects, impurities, heteroj unction roughness, and/or interaction with short- 
wavelength phonons. However, in (001) heterostructures, interaction of F valley with 
Xz valleys (i.e, X valleys lying on the quasi-momentum axis directed along the normal 
to the surface) is caused by the potential of the structure itself which leads to non- 
conservation of the quasi-momentum component kz perpendicular to the junction. As 
a result of size quantization for fairly thin (less than 50 A) AlGaAs layers, X^-valley 
states are situated below all X valleys so that an analysis of T-X^ interaction at the 
heterojunction is important for an accurate description of T-X transitions 

This type of intervalley mixing has been studied theoretically both phenomenologi- 
cally and using a tight-binding model 0, |^, and also using a pseudo-potential method 
1^. Nevertheless, the level of understanding of the processes leading to F-X^ mixing of 
electronic states is far from satisfactory. For instance, the results of Ref. p, ^ indi- 
cate that the direct Fi-Xi interaction is extremely weak and T-Xz mixing is merely 
attributable to F1-X3 interaction, whereas, according to Ref. |]^, Fi-Xi interaction is 
the determining factor (however, the parameters of the tight-binding model used in 
Ref. [0 are such that X3 band is very high in terms of energy and is in fact eliminated 



from the analysis). Studies ||10| , pJ| , in which the form of T-Xz interaction potential was 
determined by direct calculations of the matrix elements of the model heterointerface 
potential using the complete wave functions of the states also give a contradictory an- 
swer to the question of the strength of Fi-Xi interaction. According to the results of 



Ref. |]T0[, this interaction is weak whereas the results of Ref. |]TT| suggest the opposite: 
the strength of Fi-Xi mixing can be comparable to those of F1-X3 mixing. The solution 
of this problem is important for interpreting tunneling and optical experiments since 
Fi-Xi transitions are usually observed experimentally. 

In the present study F-X^ interaction (mixing) is analyzed using the method of 
envelope functions. This is a fairly explicit technique which does not have the disad- 
vantage of the tight-binding model in which the heterointerface is oversimplified, while 
calculations using the empirical pseudo-potential method are much more cumbersome 
than those using the envelope function method. Note that so far attempts to derive a 
system of equations for the envelope functions jointly describing F and X states have 
not produced satisfactory results. This is mainly because an atomically abrupt change 
in the crystal potential at the heteroj unctions produces F-X^ mixing |12[, and correct 
allowance for such abrupt changes in the potential is outside the scope of the usual 



method of Luttinger-Kohn envelope functions. In Refs. [ITU], O], in which the prob- 
lem was analyzed using the envelope function method, F-X interaction was analyzed 
using perturbation theory, but as a basis there was selected a generally over-full set 
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corresponding to the set of Kohn-Luttinger functions for F and X states. An over-full 
(and nonorthogonal) basis can, in principle, give an erroneous result. In addition, this 
approach cannot be applied directly to describe the states of the continuous spectrum 
which is important for T-Xz tunneling problem. Nevertheless, the results of the present 
study agree qualitatively with the conclusions of Ref. [Tl[ . 

The problem of adequate description of intervalley mixing of states in heterostruc- 
tures is similar to the problem of intervalley splitting of impurity states in multi-valley 



semiconductors (for example, in Si and Ge). It is known, see Ref. |13], that allowance for 
the short range part of the impurity potential ("the correction to the central cell") not 
only yields a "chemical shift" of the impurity-state energy but also lifts the valley degen- 
eracy. In Ref. [|14|, we proposed a fairly simple method of analyzing heteroj unctions with 
an atomically abrupt change in chemical composition. This method involves isolating 
the "smooth" component of the heterostructure potential and the "sharp" component 
which is only nonzero near the heterojunction. The smooth component is "processed" 
by a standard technique (the Kohn-Luttinger method) while the sharp component is 
considered as a correction to the central cell. In Ref. we only considered states near 
r point in the Brillouin zone and in the present study we develop the method further to 
describe the interaction of states near different points in k space. These results were first 



presented at the III All-Russia Conference on the Physics of Semiconductors, see [|I5 
and also |1T6[| . 



2 Formalism of the Envelope Function Method 
2.1 Formulation of the Problem 

We shall consider the electron states in (001) III-V heterostructures formed from re- 
lated lattice-matched (with the lattice constant a) semiconductors having zinc blende 
symmetry. By related structural materials we understand fairly small band offsets so 
that in the energy range of interest the conduction band states near F point can be 
described using a single-band variant and the states near the point X can be described 
using two-band (for Xi and X^ bands) variants of the envelope function method. We 
shall also assume that the energy gap between the states of interest to us in the F and 
X valleys is smaller than or of the order of the band offset. This means that we can 
consider the direct interaction of Fi, Xi, and X^ states "exactly" and interaction via 
all other bands can be taken into account using perturbation theory. Since the desired 
Hamiltonian of the equation for the envelope functions should depend explicitly on the 
number of monoatomic layers of each material forming the structure 0, [l^, we shall 
consider a structure with two symmetric heteroj unctions as the simplest case to obtain 
this dependence. For simplicity, we shall neglect spin-orbit interaction and also exter- 
nal smooth potentials. The single-electron Schrodinger equation then has the following 
form: 
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Here mo is the free electron mass, p is the momentum operator, and U = U (r) is the 
crystal potential of the heterostructure. We shall first use the following model for this 
potential (a more realistic situation will be discussed in Section 3): 

U = Ui + P{z) [U2 -Ui] = Ui + P (z) 6U, (2) 

and Ui = Ui (r) and U2 = U2 (r) are periodic (continued to all space) potentials of the 
two heterostructure materials, the z-axis is perpendicular to the hetero junction plane; 
the form-factor P{z) of the heterostructure having the heterointerfaces at 2; = and 
z = L is defined so that 

r 0, z < -d] 

P{z) = \l, d<z<L-d- (3) 
lo, z>L + d. 

The behavior of the function P{z) in the transition regions (of width 2d each) near the 
heterointerfaces may be fairly arbitrary. Here, the symmetry of the structure implies 
that P{z) = P{L — z), L > d. We shall assume that the layer of width L contains an 
integer number of monolayers: L = Afa/2, where A/" is a natural number. 



2.2 Choice of Complete Orthonormalized Set of Basis Func- 
tions 

As in Ref. |T^, the potential P{z)6U is analyzed in terms of perturbation theory and 
the complete and orthonormalized set of functions used to expand the complete wave 
function \E'(r) is constructed of Bloch functions of the base semiconductor (having the 
crystal potential Ui). In our case, the most natural basis is a mixed basis of Kohn- 
Luttinger functions for the points F and X^. In order to make this set complete and 
orthonormalized, it should be constructed as follows. We first expand \E'(r) in terms 
of the Bloch functions Mnk(r)e*^'" of the base crystal which correspond to the energy 
eigenvalues e„k, where n and k are the band index and the quasi-momentum vector, 
respectively: 

v]/ (r) = ^ 1 An' (k') e'^^'^Un'^,' d^k'. (4) 

n' 

Summation in (^) is performed over all bands and integration is performed over the 
region Aq nonequivalent k. In order not to consider two equivalent Xz points (having 
the coordinates (0,0,27r/a) and (0,0,— 27r/a) in k space), we shall not operate in the 
first Brillouin zone constructed as a Wigner-Seitz cell but we shall define a region Aq 
such that points F and, for example, q = (0,0,27r/a) are contained in this region with 
their vicinities. Following Ref. [|T^ (Sec. 7-3), we divide Aq into two subregions Ar 
and Ax containing the points F and X^ with their vicinities, where Ar U Ax = Aq and 
Ar n Ax = (see comments on this method of division in Section 2.4). Now, following 



Ref. [|T§] we use series expansions of the periodic function Unk- 
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Then can be rewritten in the following form: 

VI/ (r) = ^ / An' (k') (k') e^'^'^M^^o d3fc'+ 

+ ^ / _ A' (k' + q) c^'n' (k' + q) e'^'^e^'u^'^ d^k'. (5) 
We then define the functions 

^i'^(k')=E^n'(k')&w(k'), (6) 

•^i^^ (k') = E ^n' (k' + q) C„.'n' (k' + q) , (7) 

n' 

which will specifically comprise the envelopes of the functions of F and Xz states in k 
representation, and we shall define the region Ax such that the condition k G Ax is 
satisfied for all k + q G Ax (shift of the origin in k space). Now (^ can have the form 
of an expansion of \1' (r) in terms of a complete set of Kohn-Luttinger functions near the 
points r and X^: 

^ (r) = E / -^i'^ (kO £1 ^'k' + T^^^ (k') xS^ d^fc' = 

= E E / ^in? (k') xil d'k', (8) 

v'=r,X m' ''^v' 

where the Kohn-Luttinger functions are 

(r) _ ikr _ ikr , (r) (X) _ ikr jqr _ ikr , (X) 

Amk ~ "mO — C (p^ , Xmk ~ ^ "rrtq — C . 

Thus, the Fourier transforms of the envelope functions constructed above differ from the 
usual ones [|18] only in terms of the domains of their definition. In our case, these are 
the regions Ap (for states near the center of the Brillouin zone) and Ax (for states near 
Xz point) rather than the complete first Brillouin zone. 

For simplicity, below we shall also assume that are real. Taking the following 
orthonormalization relationship for the Bloch functions: 

/ <'k'e"''^''"Wnke*'''- dV = 5nn'S (k - k') , (9) 



we obtain the required orthonormalization relationship for the basis functions |18 



(xSO* = ^-'^""'^ - k') . (10) 
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2.3 Multi-band System of k • p Equations 

In the basis specified above the procedure for obtaining k ■ p system of equations is 
trivial, see [|l^. Using the expansion (H) in (|l|), multiplying both sides of the equation 



by (Xmk)* ^'^d integrating over all r space, we obtain the following system of equations: 

+ E e / Aii:;;;(k,k')^i";'(k')d'*'=e^i"'(k). (n) 



v'=T,X m' "'^v' 



Here e^;^) = e^o and e^^) = e„q. 



it J cell 



(m, f I p I TTi', f ') = for V 7^ f and i7 is the unit cell volume, and 
Mt:2' iKk') = f ^„ e-^('^-'^>0L^)p (.) 5f/0L^? dV. 



The matrix elements A4mm' (k, k') were analyzed in Ref. It was shown that the 
contribution of the perturbation potential can be divided into "smooth" and "sharp" 
components (the latter is exponentially small for smooth perturbations on the scale a). 
The contribution of the sharp component is a correction to that of the smooth component 
(for the case when the width of the heterostructure layers is much greater than a), and 
it can be written in the form of converging series in powers of {k^ — k'^). We will use the 
effective-mass approximation with spatially independent effective-mass parameters and 
we will only allow for the effects of abruptness of the heteroj unctions in the first order 
in terms of the parameter akz Jl^ |19|, where k^ is the characteristic quasi- momentum 
of the state. We shall consider the intervalley elements A^^m' (k, k'), f 7^ v' in greater 
detail and for Ai'Zm' (k, k'), following Ref. M, |T3, we obtain 



(k,k') 



V (kz - k-) 5U^:l, + - E ^ — — 



a a \ - 

J P'{z)e-'''^'dz - e-'('^^-'=^)^ J P'(z)e'^^Mz 6 (ky - kj,) . (12) 

-d / - 

Here we introduce the notation: V{kz) is the Fourier transform of the function P{z)\ the 
matrix element SU^l^, = {m^v \ 6U \ m',v); Kj = (47r/a)j, and j is an integer; P'{z) = 
dP{z)/dz; and ky = {kx,ky,0). We used the symmetry P{z) = P{L — z), and also the 
fact that e*^^^ = 1. Here, we give all the matrix elements ([l^ ) required subsequently. 
We assign the indices 'w', 'u' and 'v' to the conduction-band states Fi, Xi and X3, 
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respectively, and omit the valley indices since the band index in our approximation now 
uniquely defines the state. Using symmetry concepts we obtain 



M,s (k,k') 



5 (k|| - kj| 



where s = w, u, v; and the parameters dss are determined as follows 

d 



^ {S \dU cos {KjZ) \ S) f . . ^, 

dss = -2^ ^ — j P (^) sm {Kjz)(lz. 



{s I 6U cos (Kjz) I s) 

The potential also makes a contribution to the direct interaction of Xi and X3 states: 
M,,. (k, k') = i-rf,. (1 - e-^^'^-''^^"^) 6 (k|| - kj|) , 

where 



(vv') 

We shall now consider the most interesting to us intervalley matrix elements M-l^^i, 



V 7^ v': 



We shall analyze this matrix element using the same method which yielded (|T2D. We 
shall use an expansion of the periodic function UmoSUum'q as a Fourier series which gives 

(k, k') = E Cr^^^-^'^'^^V {h -K-q. + K^i) 5 (k|| - kj| + K||,) , (13) 
where K/ = (i^2i,K||/) are the vectors of the reciprocal lattice, and 

m(r)W(x) = / w^o^f/e^'^'^w^'qdV = (m, T \ 5U e'^'^^-^'^^ \ m', x) . 



As was shown in Ref. |T9||, for the region where \kx\ + \ky\ < vr/a, whose size is fairly 
large for our purposes [we obviously used this constraint in the derivation of (|12D], only 
the vectors of the reciprocal lattice with Ky; = will contribute to (pTSl): 

M^^} (k, k') = ^ (m, r I 5f/e^(^^-''^)^ I m', x) V {h -K + K, - q,) 6 (k,, - kj,) = 

j 

J2 r I 5Ue''-^^' I m', X) V (k, -k', + —j) 6 (k,, - kj,) . (14) 

j=±l,±3,±5,... ^ ^ 

We shall analyze the function V{kz — k'^ + (27r/a)j), contained in (p^ ) in greater detail: 

+ C30 

V(k,-k'^ + —j] = —- I P' {z) e-i^^-^'^+'-^^y dz = 



2 ' a 



-00 



2vr^ fc, - A;; + 



where we again used the property P{z) = P{L—z). We now expand (kz—k'^ + {27i/a)j)~^, 
exp{—i{kz — k'^)z), and exp(i(/c^ — k'^)z) as series in powers of {k^ — k'^), where the 
convergence of the first series is ensured by the property | k^ — k'^ \< 27r/a, since k G 
Ar and k' G Ax- We only retain the first terms of these expansions. This gives a 
good approximation for k^ ^ 27i/a and k^ ^ l/{2d), where k^ is the characteristic 
quasi-momentum of the state or the reciprocal characteristic length of variation of the 
envelope functions and (|^) in r representation (i.e., the envelope functions should 
vary weakly on scales of order a and on scales of the order of the widths of the interface 
regions). Now, bearing in mind that for all values of the summation index j in (p^ 
exp{—2TTijL/a) = (— 1)-'^, we obtain: 



d 



-d 



Here, there is a dependence of the effective potential on the number of monoatomic 
layers A/", which mixes the valley states assigned to different points in k space. We now 
write the matrix elements we require (we again drop the valley index): 



where 



M^. (k, k') = ^d^^ (l + (-l)-^e-^('=^-'=^)M 5 (k|| - kj| 



2 2 

dwn = - E 7r~ (w I 5f/ COS ( — jz) I u ) / P'{z) sin ( — jz) dz; (15) 
27r? \ a I J a 



i=±l,±3,±5,... 



and 



where 



(k, k') = i-rf^. f 1 - (-1)^ e''^''^-''^^'^) 5 fkii - k 



d„ 



2 d 2 

(w|5[/sin( — jz) I v) / P'(2;)cos( — jz)dz. (16) 

i=±i,±3,±5,... 27rj \ a / _7 a 

It can be seen that the abrupt potential of the heterointerfaces not only ensures Ti-X^ 
§, |T7|, |Tg, but also Ti-Xi interaction fTTJ. 

We have obtained a multi-band k ■ p system of equations and explicit expressions for 
all the elements of this system required subsequently. We shall now obtain the system 
of three equations for the strongly interacting states Fi, Xi, and X3. 
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2.4 Elimination of Remote Bands and Transition to r Space 

In the effective-mass approximation with spatially independent effective-mass parame- 
ters, the approximate unitary transformation procedure which eliminates the influence 
of remote bands in the required order of perturbation theory is performed by a standard 
method and is not given here. We obtain an integral system of equations in k rep- 
resentation. The aim of the present study is to simplify the final results as far as possible 
which may be achieved if the corresponding equations are differential. The problem of 
the accuracy of the envelope function method which occurs on transition from integral 
to differential equations and also when the remote bans are eliminated has already been 
discussed in Refs. |T^, ^ . We shall summarize the constraints imposed on the accuracy 



of the envelope function method for this case. First, we need to determine the effec- 
tive radius of the region in k space for which the system of equations for the envelope 
functions is correct (we shall call this k region the fundamental region). Whereas the 
multi-band system of equations ([Tl| ) is valid for all k G A^, and k' G A^', the unitary 
transformation of this system which eliminates the remote bands can, in principle, re- 
duce the dimensions of the fundamental k regions. This is easily understood from the 
following. The spectrum of states of a bulk semiconductor e^'^\li) near the point ko 
in a band numbered n may be represented as a series in powers of k (for degenerate 
states the spectrum is determined by diagonalizing the matrix whose elements are such 
series). The series has a finite radius of convergence Rq, which is determined by the 
strength of the k ■ p interaction with the remote bands. This radius can be estimated 
as Rq = ■moEg/{2hP), where Eg and P are the characteristic values of the interband 
energy at point kg and the interband matrix element of the momentum, respectively. 
States having quasi-momenta which do not belong to the fundamental region cannot 

be correctly taken into account in the transformed equation and should be neglected. 

fr) (X) 

We denote the corresponding radii of the fundamental k regions as R};' and i?r for 
conduction-band states near the points T and X, respectively. 

Secondly, on changing from k to r representation, integration is performed over 
regions of k space of finite dimensions which prevents us from direct obtaining differential 
equations in r space. The local approximation formula involves replacing the finite 
regions of k space by infinite ones. Since the heterostructure potential is not smooth 
and the envelope functions obtained (or their derivatives) can vary appreciably on scales 
of the order of a, this procedure does not give an exponentially small error as in the 
case of smooth perturbations, but an error which is only small in terms of power (for 
this analysis it is convenient to consider the limiting case of a mathematically abrupt 
potential and then the envelope functions or their derivatives obtained as a result of the 
local approximation may have a discontinuity). Since the perturbation-theory series used 
for the unitary transformation is a power series, it is important to avoid the inclusion 
of extra-accuracy terms. For the simple single-valley case (F states) analyzed in detail 



IS 



in Ref. [11^1 , the error of the method is of the order of (k^p / Rq )^, where l/k^ 
the characteristic length of variation of the corresponding envelope function and the 
exponent M is a measure of the smoothness of the latter (for an isolated heteroj unction or 
a fairly wide quantum well M = 3). In this case, we need to determine two characteristic 
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quasi-momenta k^f^ and fc^^-* for T and X states. In the zero-order approximation, 
which is acceptable for estimating the accuracy of the method, the F and X states do 
not interact and the constraint associated with the transition to differential equations 
is determined by the error {kP/RPf for r states (we shall assume that the layer 
width L is sufficiently large so that k^L > 1). For X states, the situation is slightly 
more complex: the system of equations for these contains both second derivatives of the 
envelope functions with respect to z and first derivatives as a result of k ■ p interaction 
between Xi and X3 bands, where the effect of this interaction may be comparable with 
the contribution of terms which are quadratic with respect to the momentum operator. 
This means that in the "worst" case, the accuracy of the local approximation for X states 
is limited by the error (k'^^'> / . All these factors allow us to consider the corrections 
which appear as a result of the abruptness of the change in the heterointerface potential 
which are small as small are the parameters ak^ and akp\ Here we also note that 
splitting the region Aq into the subregions A-p and Ax should be performed so that 

(T) (X) 

R};' and Rf' are not larger than the radii of the regions A-p and A^, respectively, and 

(T) (X) 

then Rq and Rq , not the radii Ap and Ax appear in the expressions to estimate the 
accuracy of the local approximation (this was implied above). 

As a result, the required system of differential equations for the transformed envelope 
functions Fm (r), where m = w, u, v, has the following form for the strongly interacting 
Fi, Xi, and X3 states: 

Yl C^rnm' + Vmm' (z)) Fr,i> (r) = eFm (r) . (17) 

m'=-w, u, V 

Here, T and V are (3 x 3) matrices of the kinetic and potential energies. The form of 
the matrix of the effective kinetic energy operator is known |2( 

\ 

2 

\ eI_ (p^)"^p^ + -yd d 

^ 2mu|| ^ 2mux mo ^ WxPy 

V ^ mo ^ ~^ TPa'Py 2^ + 2^ J 

where is the effective mass for F conduction band states, m^n and mi±_ are the 
longitudinal and transverse effective masses for / band (/ = u, v) at point X; and the 
bulk parameter 7, in particular, determines the magnitude of the linear photogalvanic 



effect [^. Before giving the form of the matrix of the potential energy operator, we go 
over from the function P{z) to the function Q[z) — Q{z — L) merely for the reasons of 
convenience 



19 : 



P{z)^Q (z) -Q{z- L) + poi6 (z) + 5 (z - L)) , (18) 

where 

d 

Po = J P{z) dz-d. 
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Now the matrix of the effective potential energy operator V may be expressed as the 
sum of three matrices: Vi + V2 + V3. The diagonal matrix Vi corresponds to the 



standard (bulk) effective-mass approximation |18]: 

(Vi),,, = {ess' + SUss' [6 (z) -e{z- L)]}Sss'. 

As a result of the abruptness of the heterointerface potential, the matrix V2 contains 
both intravalley contributions and contributions which mix Xi and X3 states (all the 
parameters dmm' are real): 

/d^^{6{z) + 6{z-L)) 

V2 = {6{z) + 6{z - L)) d^, {6{z) - 6{z - L)) 

\ d^^{5{z)-5{z-L)) d,, {6{z) + 6{z - L)) 

Here, the parameters dss are related to dss as follows: 

dss = dss + SUssPo- 

Finally, the matrix V3 contains contributions which mix F and Xz states (also as a result 
of the abruptness of the heterointerface potential). We give the non-zero elements of 
V3: 

' {z) + (-l f6(z-L)], (19) 

(20) 



= (V,)^^ = {6{z) + {-lf6{z - L)) , 

In our approximation, the expression linking the envelope functions and the complete 
wave function has the usual form (we again drop the valley index): 



E 



m=w, u, V 



+ E 



imo {em - en) 



For an arbitrary number of heteroj unctions at z = Zj numbered by the index j, the 
nonzero elements of V3 may be written in the following form: 



(Vs. 



(Vs. 



(Vs. 



z — z 



3I1 



where q is the distance in k space between the centers of the valleys under consideration 
(in our case, q = (0, 0, 27r/a)). Since the parameters d^l and d^l depend not only on 
the heteroj unction materials but also on the microscopic structure of the interface, these 
must generally be determined for each interface separately. 
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3 Discussion of Results and Conclusions 



Using k ■ p formalism, we have constructed a generalization of the method of envelope 
functions suitable to describe the interaction of F and Xz states in (001) III-V nanos- 
tructures formed from related lattice-matched semiconductors. In the derived system of 
equations ([TtD , mixing of the states of different valleys is determined by the heterosur- 
face effective potentials (p!9| ) and (pOD, similar to those introduced phenomenologically 
in Ref. P]. The system (|l^) contains information on the number Af of monoatomic 
layers of the structure. An oscillatory dependence of the effective intervalley mixing 
strength on Af was obtained, in particular, in Ref. [0], and was also introduced from 
symmetry concepts (in terms of the boundary conditions for the envelope functions) in 



Ref. [|T^. However, in addition to this, there are some difference between the results 
of the present study and the results of other authors. The most important of these is 
the appearance of direct interaction between Fi and Xi states whose strength is deter- 
mined by the parameter d„^. This was predicted recently in Refs. fll], |15|, [l6l. The 
strength of Fi-Xi interaction, i.e., the value of the constant c/^u; depends strongly on 
the structure of the heterointerface on atomic scales. In simplified models, such as the 
simplest variants of the tight-binding method, this interaction may be absent. It is clear 
from (|15D that in the hypothetical case of mathematically abrupt heteroj unctions when 
P{z) = Q{z) — Q{z — L), we in fact find d^^ = 0. The conclusion reached in Ref. |]TU| that 
dwu d^v was specifically a consequence of the selected heteroj unction model for which 
the simplest case is a mathematically abrupt jump (see also Ref. [|TT|). In general, there 
is no basis for assuming that d^u differs substantially from d^^, and since Xi band is 
lower in energy than X3 one, the existence of Fi-Xi interaction may be very significant 
for interpreting experiments. 

In Ref. we considered a more complex model of the potential of an ideal het- 
erostructure and, specifically, took into account a periodic coordinate dependence of the 
heterostructure form-factor in the (001) plane. We showed that this complication merely 
leads to renormalization of some parameters obtained using (Q). However, whereas it 
was found using the simple form-factor model that the mixing strength of the heavy and 
light holes at the center of the 2D Brillouin zone was higher for abrupt heteroj unctions 
than for structures with a continuously varying composition, this conclusion was not 
obtained using the more complex form-factor model. We could merely conclude that the 
strength of the light hole-heavy hole mixing depended on the structure of the transition 
region of the heterointerfaces. Having made a similar analysis for our case, we can easily 
show that F-X mixing strength in fact depends strongly on how abruptly (on scales 
of the order a) the chemical composition at the heteroj unction varies. This conclusion 
which follows directly from expressions (p^ ) and (|T6]) for the parameters rf^u and d^v is 
quite natural since transitions with such a large (of the order of the dimensions of the 
Brillouin zone) change in the quasi-momentum kz can only be achieved by electron scat- 
tering at an atomically abrupt heterointerface (in formal terms the Fourier components 
of the heterointerface potential with the wave vectors close to kz = 2tt/ a are responsible 
for F-X transitions). 
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